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ABSTRACT 


Methods  are  described  for  obtaining  similarity-law  characteriza- 
tions for  arbitrary  roughnesses  indirectly  from  measurements  of  pipe- 
flow  head  loss,  flat-plate  drag,  and  rotating-disk  torque.  The  rough- 
ness characterizations  are  used  to  make  viscous-drag  predictions  by 
means  of  boundary- layer  calculations.  A simplified  drag-prediction 
method  using  form  factors  and  equivalent  rough  flat  plates  is 
described  in  detail. 

ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  funded  by  the  Naval  Sea  Systems  Command 
(PMS  393),  Task  Area  S0411001,  Element  64561N. 

INTRODUCTION 

The  prediction  of  the  increase  in  drag  due  to  a particular  roughness 

on  ships  does  not  seem  to  be  a straightforward  procedure,  especially  for  a 

roughness  with  an  arbitrary  geometrical  configuration.  One  paper  advocates 
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just  full-scale  ship  tests  while  another  considers  only  deviations  from 
smooth  conditions  for  a towed  friction  plane  with  various  roughnesses  without 
further  analysis  for  full-scale  conditions.  The  role  of  the  similarity-law 
roughness  characterization  does  not  seem  to  be  understood  in  extrapolating 
results  obtained  under  model-scale  conditions.  The  prediction  of  the 
viscous  drag  of  rough  bodies  by  any  boundary-layer  analysis,  even  at  the 
equation-of-motion  level,  requires  as  an  input  a similarity-law  characterization 
for  the  particular  roughness  configuration  being  considered. 

At  high  Reynolds  numbers,  the  primary  effect  of  roughness  is  in 
increasing  turbulent  skin  friction  which  is  the  case  to  be  considered  here. 

^Johnsen,  S.,  "Drag  Reducing  Coatings  for  Ship  Hulls"  presented  to  8th 
2Skandinaviska  Lackteknikers  Forbund  Congress,  29  Sep  - 1 Oct  1976. 

West,  E.E.,  "The  Effect  of  Surface  Prepaiation  and  Repainting  Procedures 
on  the  Frictional  Resistance  of  Old  Ship  Bottom  Plates  as  Predicted  from 
NSRDC  Friction  Plane  Model  4125,"  Naval  Ship  R8D  Center  Report  4084 
(May  1973) . 


Since  the  physics  of  turbulent  flow  and  turbulent  skin  friction  is  not 
sufficiently  developed  to  provide  theoretical  predictions  for  even 
smooth  surfaces,  let  alone  rough  surfaces,  the  velocity  similarity  laws 
have  proved  immensely  valuable  in  supplying  engineering  solutions  on  a 
semi -empirical  basis. 

The  velocity  similarity  laws  have  proved  to  be  very  successful  in 

correlating  turbulent  shear  flows,  internal  and  external,  such  as  pipe  flow 
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and  boundary-layer  flow,  and  for  smooth  and  for  rough  surfaces.  ’ 

For  rough  surfaces,  the  velocity  similarity  laws  lead  to  roughness 
characterizations  which  have  been  verified  in  principle  by  the  classical 
experiments  by  Nikuradse^  on  dense  sand-grain  roughness.  Here,  different 
size  sand  grains  in  pipes  of  different  diameters  gave  friction  losses  and 
velocity  profiles  which  are  closely  correlated  by  a similarity-law 
characterization.  Hama^  has  shown  that  a roughness  characterization  is  the 
same  for  both  pipe  flow  and  boundary- layer  flow. 

An  inherent  problem  with  irregular  roughness  is  the  difficulty  of 
actually  manufacturing  a reduced-scale  roughness  if  desired  to  suit  reduced- 
scale  laboratory  experiments.  The  most  feasible  procedure  is  to  treat  an 
irregular  roughness  as  a full-scale  roughness  per  se  even  in  laboratory-scale 
experiments . 

A direct  roughness  similarity-law  characterization  requires  the 
measurement  of  a velocity  profile  correlated  with  a local  wall  shear 
stress.  This  is  usually  an  arduous  and  difficult  undertaking,  especfially 
when  the  shear  layer  is  thin.  It  is  much  easier  to  obtain  a roughness 
similarity-law  characterization  indirectly  from  more  convenient  overall 
measurements  such  as  the  average  velocity  and  pressure  drop  in  pipe  flow, 
the  drag  of  a towed  flat  plate  and  the  torque  of  a rotating  disk.  The 
indirect  characterizations,  however,  require  appropriate  analytical 
relations  derived  from  the  similarity  laws  themselves. 


^Schlichting,  H.,  "Boundary-Layer  Theory,"  6th  Ed.,  McGraw-Hill,  New  York, 
.1968. 

gWhite,  F.M.,  "Viscous  Fluid  Flow,"  McGraw-Hill,  N.Y.  (1974). 

^Hama,  F.R.,  "Boundary-Layer  Characteristics  for  Smooth  and  Rough  Surfaces," 
Transactions  of  Society  of  Naval  Architects  and  Marine  Engineers,  Vol . 62, 
pp.  333-358  (1954). 
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The  purpose  of  this  paper  is  to  present  and  assemble  explicit 
procedures  for  obtaining  roughness  similarity-law  characterizations  in- 
directly from  overall  measurements  of  pipe  flows,  of  towed  flat  plates, 
and  of  rotating  disks  based  on  appropriate  analytical  relations.  The 
presentation  also  includes  methods  of  utilizing  the  roughness  characterizations 
to  make  predictions  of  the  viscous  drag  of  undeiivater  bodies  from  various 
boundary- layer  analyses.  Emphasis  is  given  to  a simplified  method  based  on 
form  factors  and  equivalent  flat-plate  drag  charts  for  particular  roughnesses. 
Simple  methods  are  presented  for  preparing  these  flat-plate  drag  charts  from 
similarity-law  characterizations.  Numerical  examples  are  given  to 
illustrate  the  procedures. 

VELOCITY  SIMILARITY  LAWS 

The  present  state  of  development  of  the  velocity  similarity  laws 
for  turbulent  shear  flows  past  irregular  rough  surfaces  are  first  reviewed.^ 

The  classical  velocity  similarity  laws  developed  by  F’randtl  and  by 
von  Karman  for  flow  through  smooth  pipes  led  to  a logarithmic  formula 
which  successfully  correlated  existing  frictional  data.  The  similarity 
laws  were  extended  to  rough  pipes  by  Nikuradse  who  successfully  correlated 
the  frictional  data  of  dense  sand-grain  roughnesses  of  different  sizes 
in  various  diameter  pipes.  A variety  of  other  roughness  configurations 
have  also  been  investigated  since  then.  In  time,  the  v'elocity  similarity 
laws  have  been  extended  to  other  shear  flows  like  boundary-layer  flows. 

In  general  for  turbulent  shear  flows  the  two  laws  which  provide 
similarity  to  the  mean  velocity  profile  by  relating  it  to  the  wall  shear 
stress  are: 

1.  The  inner  law  or  the  law  of  the  wall  which  applies  to  the  flow 
immediately  adjacent  to  the  solid  boundary. 

2.  The  outer  law  or  the  velocity  defect  law  which  applies  to  the 
remaining  outer  region  of  the  shear  flow. 

As  originally  analyzed  by  Millikan,  consideration  of  an  overlapping  of  the 
two  similarity  laws  in  the  shear  flow  results  in  a logarithmic  functional 

"’’^e  Erictional  Resistance  and  '^urbtilent  Boundary  Layer 
of  Rough  Surfaces,"  Journal  of  Ship  Research,  Vol . 2,  No.  3,  pp.  52-74 
(Dec  1958). 
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form  for  both  similarity  laws  within  the  common  region  of  overlap.  Outside 
the  overlapping  region,  a functional  form  is  provided  by  the  law  of  the 

7 

wake  and  a wake  modification  function. 

Inner  Law  or  Law  of  the  Wall 

A similarity  law  may  be  developed  for  the  mean  velocity  u of  the 
turbulent  flow  in  the  region  close  to  a rough  wall  with  an  irregular 
configuration  by  considering  the  roughness  of  the  wall  defined  by  any 
number  of  length  parameters  k,  k^,  k^,  ...as  an  addition  to  the  usual 
parameters  of  wall  shear  stress  density  p and  kinematic  viscosity  of 
the  fluid  *>,  and  normal  distance  from  the  wall  y or 


By  dimensional  analysis,  the  variables  may  be  grouped  into  the  following 
nondimensional  ratios: 


or 


u 

U-c 


where 


(2) 


(3) 


(4) 


U.J  is  the  shear  velocity^  U-j  and 


D c k-  ^ 


(5) 


(6) 


^Granville,  P.S.,  "Similarity-Law  Entrainment  Method  for  Two-Dimensional 
Turbulent  Boundary  Layers  in  Pressure  Gradients,"  David  Taylor  Naval 
Ship  R&D  Center  Report  4657  (Dec  1975). 
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are  dimensionless  ratios  representing  the  shape  of  an  irregular  roughness 
configuration . 

Outer  Law  or  Velocity  Defect  Law 

In  general  for  turbulent  shear  flows,  the  velocity  defect  U-u 
has  been  experimentally  found  to  be  directly  independent  of  viscosity, 
except  close  to  the  wall  and  a function  only  of  wall  shear  stress  , 
density  of  fluid  |»  , and  distance  "nwardi-y  or 


U-U  ^ f 


(7) 


Here  U is  the  velocity  at  the  center  of  a pipe  or  the  velocity  at  the 
outer  edge  of  a boundary  layer  of  thickness  S . 

By  dimensional  analysis 


0- 


a 


(8) 


represents  the  outer  law  or  velocity  defect  law  which  has  been 
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experimentally  found  to  be  unaffected  by  a rough  wall.  ' 

Logarithmic  Velocity  Law  and  Similarity-Law  Roughness  Characterization 
Both  the  inner  and  the  outer  similarity  laws  may  be  considered  to 
be  effective  in  a common  overlapping  region  of  a shear  flow.  Equating  the 
y-derivatives  of  the  inner  and  the  outer  laws  results  in  a logarithmic 
relation  for  both  laws  in  the  overlapping  region.^ 

For  the  inner  law,  the  logarithmic  velocity  law  becomes 


(9) 


where 


A = constant 
• 'l,0  * 

Bj  = constant  value  of  B^  for  smooth  walls 


(10) 


and 


(11) 
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which  may  be  termca  The  ■ 

For  smooth  walls  AB  Tti  • • . • 

from  the  overlapping  o'  Thv  • • t 
character! tat  ion  may  hi  jr:^.  n n 
and  is  a key  element  ir.  ri.r 
hydrodynamic  bodies  due  tc  r >i»t  • - 
AB  is  negative. 

An  alternate  t•xpre^s.or 

^ = A.  u ^ . 

' k 

where 

B,  = b.  - Am-* 

For  the  case  of  fulU  devr  • .. 

value  and 

Ab  ■ Br,o  - - A .r  ■* 

is  then  a straight  line  ■ r ♦. 
of  AB  against  In  U*. 

For  the  outer  law,  thr  • • 


U - u 


which  is  the  sane  a«  'hat  tor 
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A smooth  surface  gives 
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At  the  same  y then 


(16) 


(17) 
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The  direct  velocity  method  was  used  by  Nikuradse  in  his  classical 
flow  measurements  for  dense  sand-grain  roughness  in  a pipe  flow.  Since 
such  velocity  measurements  may  be  difficult  and  arduous  to  perform,  it 
is  more  convenient  to  resort  to  indirect  methods  requiring  just  overall 
measurements,  such  as  the  average  velocity  and  pressure  drop  for  pipe 
flow,  the  total  drag  of  a flat  plate,  and  the  torque  of  a rotating  disk. 

Pipe  Flow 

In  fully-developed  pipe  flow,  it  is  more  convenient  to  determine  the 
average  velocity  from  the  rate  of  mass  flow  than  to  perform  a velocity 
survey.  The  similarity  laws  may  then  be  used  to  relate  the  AB-characterization 
to  the  average  velocity. 

In  general,  for  straight  circular  pipes  with  fully-developed  turbulent 
flow  (no  entrance  effects),  the  friction  loss  or  pressure  drop  represented 
by  wall  shear  stress  may  be  related  to  an  average  velocity  V of  the  pipe 
flow,  pipe  diameter  D,  fluid  properties  of  density  p and  kinematic  viscosity 
V,  and  a geometrical  description  of  the  irregular  rough  wall  k,  k^,  k^,..., 
which  symbolically  becomes 


(19) 


A dimensional  analysis  results  in  a friction  factor  f as  a function 
of  two  dimensionless  ratios,  a Reynolds  number  R„  and  a relative  roughness 
D/k  for  a particular  roughness  configuration  R or 
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(20) 


where 

the  Fanning  friction  factor  f is  given  by 

- 

'ifV‘ 

and  the  pipe  Reynolds  number  by 

The  Darcy- Wei sbach  friction  factor  A may  also  be  used  where 


I - 4^ 


(23) 


The  AB-characteri zation  may  be  related  to  overall  factors  as 
follows.  For  a circular  pipe,  the  average  velocity  V is  given  by 

»-n  ^ ^ <^y 

where  6=D/2,  pipe  radius 

which  in  inner  law  variables  becomes 


(24) 


Application  of  both  velocity  similarity  laws  and  neglect  of  the  very 
thin  laminar  and  buffer  sublayers  results  in 


V 


B,  * 


^ AB  - 


tl 

30 


(25) 


(26) 
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Since  by  definition 


}L  - fX 

Ut  " 1 ^ 


xVf  (2, 


then  for  rough  surfaces 


V?  i YT 


For  smooth  surfaces  AB=0  and  then 


(■J?)  *1/2  '"Vf 

^ rf  o 

For  the  same  value  of  Yr  Kt 


:h-)  4 l«Vf  B„o  &.  - T ^ 


AB  = 


and  the  corresponding  k is  given  by 


which  follows  from  definitions. 


The  procedure  for  determining  AB[k*,R]  from  pipe  flow  is  then 
as  follows. 

1.  For  roughness  configuration  R inside  a pipe  of  diameter  [)  plot 
the  measured  data  reduced  to  dimensionless  ratios  f and  R^  in 
coordinates  versus  log  R^^.  Then  plot  the  smooth  friction  line 

in  the  same  coordinate  system. 


o 


2.  At  specific  values  of  abscissa  Yf  Rq  determine  AB  from 
Equation  (31)  and  k*  from  Equation  (32) . 

Flat  Plate 

For  the  case  of  a towed  flat  plate  with  an  irregular  rough  surface, 
it  is  certainly  more  convenient  to  measure  the  drag  and  towing  velocity 
than  to  measure  the  local  velocity  profile  in  a thin  boundary  layer  and 
the  accompanying  local  wall  shear  stress.  Indirectly  the  similarity 
laws  may  then  be  used  to  determine  the ^B-characterization  as  a function 
of  roughness  Reynolds  number  k*  from  the  total  drag  and  towing  velocity. 

In  general,  the  drag  D of  a rough  flat  plate  towed  parallel  to  its 
motion  depends  on  forward  velocity  U,,,,  surface  area  S,  length  L,  fluid 
properties  of  density  p and  kinematic  viscosity  v and  a geometrical 
description  of  its  irregular  rough  wall  k,  k^,  k^,...  which  symbolically 
becomes 


o = i \y->*  w,  (33) 

A dimensional  analysis  results  in  a drag  coefficient  Cp  as  a function 
of  two  dimensionless  ratios,  a Reynolds  number  R,  and  a relative  roughness 
L/k  for  a particular  roughness  configuation  R or 


(54) 


where  drag  coefficient 


(35) 
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and  Reynolds  number 


(36) 


The  AB-characterization,  may  be  obtained  indirectly  from  drag 
measurements  of  a rough  surface  by  considering  deviations  from  a 
smooth  friction  line  as  follows. 

An  application  of  the  similarity  laws  results  in  a logarithmic 
formula  for  the  drag  of  rough  flat  plates 


In  Rj. 


■{I  \ 

■ 


A 


A 


\vj  2 D, 


(37) 


Here  is  a boundary- layer  constant  defined  in  Reference  6 and  (AB)^ 
is  the  value  of  AB  for  the  boundary  layer  at  the  end  of  the  flat  plate. 
For  a smooth  wall  AB  = 0.  Then 

\ o p „ ' . 6a  B.,o  1 t\ 

. 1 

Consider  a plot  of  “nr”  against  log  Rj^Cp  for  both  rough  and  smooth 

plates.  Then  at  the  saml  values  of  R, C 

L«  r 

(<^8),  -Ifll 


The  corresponding  at  the  end  of  the  flat  plate  is  obtained  as  follows: 
From  definitions 


(40) 


(41) 
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The  procedure  for  determining ^B[k  R]  is  then  as  follows: 

9 

1.  For  roughness  configuration  R on  a flat  plate  of  length  L plot 
the  measured  data  (Cp,  Rj^)  in  coordinates^:  against  log 

the  smooth  friction  line  in  the  same  coorAnate  system. 

2.  At  specific  values  of  abscissa  Rl*“F'  tJctfirmine  AB  from  Equation 
(39)  and  k from  Equations  (40)  and  (41) . 

Rotating  Disk 

Although  the  boundary- layer  flow  on  a rotating  disk  is  three  dimensional 
in  contrast  to  the  two-dimensional  flow  on  a flat  plate  or  through  a pipe, 
the  rotating  disk  provides  a convenient  method  of  characterizing  roughness 
effects  particularly  at  high  speeds.  Determining  the  AB-characterization 
directly  by  a velocity  survey  of  the  boundary  layer  of  a rotating  disk  is 
physically  feasible  but  definitely  more  arduous  than  determining  AB 
indirectly  from  overall  measurements. 

In  general  the  torque  or  moment  M of  one  side  of  a rotating  disk  in  an 
unbounded  fluid  medium  depends  on  the  angular  velocity  w,  the  radius  R, 
fluid  properties  of  density  p and  kinematic  viscosity  v,  and  the  geometry 
of  the  irregular  rough  wall  k,  k^,  k2...  or  symbolically 

NA  = 1 (42) 

A dimensional  analysis  results  in  a moment  coefficient  C as  a 

m 

function  of  two  dimensionless  ratios,  a Reynolds  number  R_  and  a relative 
roughness  R/k  for  a particular  roughness  configuration  R or 


where  moment  coefficient 


and  Reynolds  number 


4 M 


(44) 


(45) 


12 


The /^B-characteri zation  may  be  obtained  indirectly  from  overall 
measurements  as  follows.  The  application  of  the  similarity  laws  results 
in  a logarithmic  formula  for  rough  surfaces 

where  (AB)^  and  are  the  values  ofAB  and  (AB)  ' at  the  edge  of  the 

disk,  and 


(4  6) 


(ab) 


a 


(47) 


For  a smooth  surface  AB  = 0 and  (AB)  ' = 0. 
Then 


yiOTT 


(4  8) 


Consider  a plot  of-^srt"  against  log  for  both  rough  and  smooth 


surfaces.  Then  for  the  same  value  of  RpiCw, 

K 


5* 


(49) 


The  corresponding  value  of  k at  the  edge  of  the  disk  is  obtained 
as  follows  ® 

From  definitions 


_ Ux 


(jO  R. 


R» 


(SO) 


8. 


Granville,  P.S.,  "The  Torque  and  I'urbulent  Boundary  uayer  of  Rotating 
Disks  with  Smooth  and  Rough  surfaces,  and  in  Drag-Reducing  Polymer 
Solutions,"  .loumal  of  .''hip  T’esearch,  Vol . 1'!’,  No.  1,  pp.  l^l-l'’?  (Hec  1"7S) 
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and  the  corresponding  local  friction  at  the  edge  of  the  disk  is 


The  procedure  for  determining  [k*,R]  is  then  as  follows. 

1 , For  roughness  configuration  R on  a disk  of  radius  R plot  the 

t 

measured  data  (C^,  R^^)  in  coordinates  against  log  R^^  r^- 

Plot  the  SI  oth  friction  line  in  the  sam^  coordinates. 

2.  At  specific  values  of  abscissa  Rj^  V^m'  determine  AB  from 
Equation  (49)  and  k*  from  Equations  (50)  and  (51).  The  effect  of  (AB) 
may  be  obtained  by  reiteration.  Application  to  enclosed  rotating  disks 
is  shown  in  the  Appendix. 

A check  on  the  validity  of  the  similarity  laws  for  rotating  disks 

is  given  by  the  excellent  agreement  of  the  torque  of  a rota.ing  disk 

covered  with  dense  sand  grains  with  similarity-law  predictions  for  the 
9 

fully  rough  regime  . 


(SI) 


1 


APPLICATION  OF  SIMILARITY-LAW  ROUGHNESS  CHARACTERIZATIONS 
TO  PREDICTING  VISCOUS  DRAG 

Boundary-Layer  Calculations 

Since  the  viscous  drag  of  bodies  may  be  determined  from  the  boundary- 
layer  development,^^  the  increased  drag  due  to  roughness  may  also  be 
determined  by  applying  the  roughness  characterization  to  a turbulent 
boundary- layer  prediction  method.  In  general  there  are  now  two  principal 
turbulent  boundary- layer  methods: 

1.  "Differential"  methods  for  solving  the  boundary- layer  equations 
of  motion  directly  by  using  a turbulent  shear  stress  model  such  as  that  of 
eddy  viscosity  or  mixing  length. 

9 

Dorfman,  L.A.,  "Hydrodynamic  Resistance  and  the  Heat  Loss  of  Rotating 
.^Solids,"  Oliver  and  Boyd,  Edinburgh,  1963. 

^Granville,  P.S.,  "The  Calculation  of  the  Viscous  Drag  of  Bodies  of 
Revolution,"  David  Taylor  Model  Basin  Report  849  (Jul  1963). 
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2.  "Integral"  methods  for  solving  the  integrated  boundary- layer 

equations  of  motion  by  using  relations  for  velocity  profile  and  local 

skin  friction.  Tc  date,  no  "differential"  method  has  been  published 

which  includes  roughness  effects  although  Cebeci  and  Smith^^  allude  to  a 

possible  method.  The  effect  of  roughness  has  been  included  in  an 

"Integral"  method  which  extends  the  two-dimensional  empirical  Head 

12 

entrainment  method  by  altering  the  local  skin  friction.  The  effect  of 

roughness  may  also  be  readily  included  in  a similarity-law  entrainment 

7 13 

method  for  two-dimensional  and  axisymmetric  boundary  layers.  ’ 

Form  Factors  for  Rough  Surfaces 

For  streamlined  bodies  for  which  no  appreciable  flow  separation 

14 

occurs  a simplified  method  has  been  developed  for  separating  the  effect 
of  shape  and  of  roughness  on  viscous  drag.  A form  factor  is  defined  which 
accommodates  the  effect  of  shape  while  an  equivalent  flat  plate  is  defined 
which  takes  care  of  roughness  effects. 

In  general  the  drag  of  a body  depends  on  the  forward  speed  Ll^,  the 
surface  area  S,  the  density  p and  kinematic  viscosity  ^ of  the  fluid,  the 
forward  length  L,  and  other  lengths  defining  the  shape  , L^,...  and  the 
geometry  of  the  irregular  roughness  k,  k^,  k^...  or 


A dimensional  analysis  results  in 

c.  = f 

where  drag  coefficient  is  given  by 


Cebeci,  T.  and  A.M.O.  Smith,  "Analysis  of  Turbulent  Boundary  Layers," 
^_Academic  Press,  New  York,  1974, 

“^Dvorak,  F.A.,  "Calculation  of  Turbulent  Boundary  Layers  on  Rough  Surfaces 
in  Pressure  Gradient,"  AIAA  Journal,  Vol . 7,  No.  9,  pp.  1752-1759 
(Sep  1969). 

oranville,  P.S.,  "Similarity-Law  Entrainment  Method  for  Thick  Axisymmetric 
Turbulent  Boundary  Layers  in  Pressure  Gradients,"  David  Taylor  Naval  Ship 
R6D  Center  Report  4525  (Dec  1975) . 
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(54) 


Cv  - 


the  Reynolds  number  by 


Dy 

ipUiS 


G.- 


UoL 


the  shape  of  the  body  by  length  ratios 

? * JL  L, 

^ * T“  ) -r-  ) • • • (56) 

U,  L J, 

and  the  shape  of  the  irregular  roughness  by  length  ratios 

A form  factor  f is  defined  so  that 

where  Cp  is  the  equivalent  flat-plate  drag  coefficient.  An  equivalent 
flat  plate  is  defined  as  a flat  plate  with  the  same  surface  area  and 
length  as  the  body. 

14 

By  a simplified  boundary-layer  analysis  it  may  be  shown  that  a form 
factor  f may  be  determined  approximately  from  the  pressure  along  the  body 

represented  by  velocity  ratio  — . Here  U is  the  velocity  outside  the  boundary 

1 ^*0 
layer. 

14 

For  two-dimensional  foils  at  zero  angle  of  attack 


I +•  " 


,1 

J'oJ 

I,  »e»-4  JIl") 


Granville,  P.S.,  "A  Prediction  Method  for  the  Viscous  Drag  of  Ships 
and  Underwater  Bodies  with  Surface  Roughness  and/or  Drag-Reducing 
Polymer  Solutions,"  David  Taylor  Naval  Ship  R6D  Center  Report  SPD-797-01 
(Oct  1977). 
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Here  H is  a two-dimensional  shape  parameter 
X is  a chord  distance  from  the  nose 
q is  a wake  factor 

e is  a subscript  representing  conditions  at  the  tail  end  and 
oC  is  an  angle  of  the  tangent  of  the  contour  relative  to  the  chord 

14 

For  bodies  of  revolution  in  axisymmetric  flow 


I 


flaia 

IB 


iec  d 


tec  d d 


lx] 


(60) 


Here 


and 


is  an  axis>Tranetric  shape  parameter 

is  a radius  of  cross-section  of  the  body  of  revolution 
is  an  axial  distance  from  the  nose 

is  an  angle  of  the  tangent  of  the  meridian  contour  relative  to 
the  axis 

12 

For  ship  hulls  two  equivalent  bodies  of  revolution  are  used  in 
Equation  (60)  to  give 


h 

r 

X 

oL 


- 


£‘t 

IB 

tec  *1 

Here  r is  a radius  of  a cross-section  of  an  equivalent  body  of 


revolution  based  on  perimeter. 


is  a velocity  ratio  of  an 


(61) 


equivalent  body  of  revolution  based  on  cross-sectional  area  defined  by  r^. 
Also 


r 


P 


Tl 


(62) 
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and 


where 

P is  the  local  perimeter  of  a ship's  underwater  hull. 

A is  a cross-sectional  area  of  a ship  underwater  hull  at  an  axial  station. 

A constant  value  of  m sucn  as  m = n = 1.4  is  recommendea  although  other 

values  of  H could  accommodate  any  changes  in  form  factor  f due  to  viscous 
effects.  For  two-dimf\^sional  foils  q = 1 and  for  bodies  of  revolution 
q = 7. 

FLAT-PLATE  DRAG  DIAGRAMS 

Flat-Plate  Drag  Diagrams  from  Roughness  Characterizations 

For  a given  irregular  roughness  configuration  R the  roughness 

characterization  is  a function  of  a single  dimensionless  ratio,  Reynolds 

* 

number  k . The  requirement  is  to  convert  AB  to  a flat -pi ate  drag 
coefficient  Cp  as  a function  of  two  dimensionless  ratios,  Reynolds  number 
Rj^  and  relative  roughness  L/k.  A convenient  method  to  determine  each 
point  Cpflog^P  Rj^,  L/k]  is  what  may  be  termed  the  method  of  the  interaction 
of  two  loci:  one  locus  to  satisfy  AB  and  another  to  satisfy  k . These 
may  be  performed  graphically.  The  first  locus  to  satisfy^B  depends  on  a 
smooth  friction  line^^  such  as 


(C  \ - O OTTle  . ^ 

-..«)*  R. 

which  agrees  with  the  well-known  Schoenherr  line  at  high  Reynolds  numbers 

but  is  more  accurate  at  low  Reynolds  number.  From  Equation  (37)  for  the 

same  values  of  C_  and  in  common  logarithms 
r 

^^Granville,  P.S.,  "The  Drag  and  Turbulent  Boundary  Layer  of  Flat  Plates 
at  Low  Reynolds  Numbers,"  Journal  of  Ship  Research,  Vol . 21,  No.  1, 
pp.  30-39  (Mar  1977). 
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(6S) 


(66) 


Hence,  for  constant  the  locus  for  the  friction  line  with  roughness 
is  offset  by  a constant  amount  in  the  log^^  Rj^  direction. 

The  second  locus  which  satisfies  k*  is  obtained  by  taking  the 
common  logarithm  of  Equation  (40)  and  by  utilizing  Equation  (42)  so  that 

^ ^ * '°3“’  T 

Each  value  of  relative  roughness  L/k  produces  another  locus. 

The  intersection  of  the  two  loci  produces  a point  (Cp,  log^Q 

for  a particular  L/k  and  R.  Other  points  are  determined  to  define  a 

friction  line  C [log  R,  , L/k]  for  roughness  configuration  R. 

10  L 

Roughness  Friction  Lines  for  Different  Flat-Plate  Lengths 

If  a friction  line  for  a particular  roughness  R is  known  for  one 
length,  it  is  often  necessary  to  obtain  friction  lines  for  other  flat-plate 
lengths.  To  find  one  way  would  be  to  determine  the AB-correlation  and 

r 

use  the  method  of  two  loci  previously  described.  A more  direct  procedure  is 
now  described  which  also  involves  two  loci. 

Each  point  on  a given  roughness  friction  line  for  a plate  of  length  L, 

* 

represents  a value  of  ^B  and  k . The  object  is  to  obtain  a point  on  a new 

roughness  friction  line  for  a different  length  plate  with  the  same  values 
* ^ 
ofi^B  and  k . The  locus  of  a line  with  the  same^B  is  given  by  Equation  (37) 

1 0(^,6  = - 

which  shows  that  the  locus  is  offset  a fixed  distance  from  the  smooth 
friction  line  in  the  log^Q  Rj^  direction. 

The  other  locus  for  the  same  value  of  k but  with  a different  length 
is  obtained  as  follows.  Consider  a combined  factor  L*  where  by 
definitions 
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For  the  same  values  of  C^,  and  sane  * " 


Here  and  are  defineci  ‘ct  . j- 

The  procedure  is  to  draw  a reterer.^r  . 

to  satisfy  Equation  (68).  Throuft  p<  * 

a locus  is  drawn  which  is  offset  a . • 

direction.  The  locus  for  the  sar-e  • j 

1 

constant  distance  given  by  loEjj  |— 
constant  L and  that  for  constant  ► • r 
line  corresponding  to  leneth  1.,  a<^  vt  „• 

MIMRl  M 

To  illustrate  the  procedure'^  de«.<  r:’r 
examined.  The  first  case  to  be  cons; jr-r 
rotating  disk  with  a random  rouchnes^ 
about  0.005  inch.^^  The  C , R«  lat*  as. 

TT  • 

against  log  R„*/^  in  Figure  1 for  fotr  .,  • 

values  of  Rj^'^C^,  the  roughness  charatte-  ••  J 
Equation  (40)  with  (AB) ^ considered  rerr  -•  • 
is  obtained  from  Equations  (50;  and  'i)  an. 

^^Amfilokhiev,  W.B.  and  A.M.  rerguson,  . a < 

Surface  Roughness  in  Dilute  Polymer  x',  • • 
Department  of  Naval  Architecture  fxpenarr’  . 
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Figure  2.  The  slope  of  the  correlating  line  gives  the  values  of  (AB) ' 
which  is  then  used  to  recalculate  the  values  of  AB  as  given  by  Equation 
(49).  The  reiterated  final  results  are  also  plotted  in  Figure  2. 

The  preparation  of  a flat-plate  roughness  diagram,  Cp,  R^,  is  shown 
in  Figure  3 for  L = 300  feet  and  L = 600  feet.  The  curve  for  the  smooth 
surface  is  plotted  from  values  calculated  from  Equation  (64) . Selected 
points  (LP,  k*)  of  the  roughness  correlation  of  Figure  2 are  treated  as 
follows.  In  accordance  with  Equation  (65)  a locus  of  constant  A B is  drawn 

K g 

offset  a value  of  - — - from  the  smooth  curve  (A  = 2.5).  A locus  of  constant 
k*  for  L = 300  feet  and  L = 600  feet  are  plotted  to  satisfy  Equation  (60) 
by  calculating  as  a function  of  Cp.  The  intersections  give  points  on 
the  lines  for  L = 300  feet  and  L = 600  feet  respectively.  The  procedure 
is  repeated  for  other  values  ofAB  and  k*.  The  resulting  curves  are  plotted 
in  Figure  3 within  the  limits  of  the  experimental  data.  The  largest  value  of 
Reynolds  number  Rj^  for  L = 300  feet  is  R^^  = 1.8  x 10  which  corresponds  to 
a speed  of  about  43  knots  which  is  higher  than  the  operating  velocity  of 
most  vessels.  This  shows  that  the  rotating  disk  easily  provides  results 
applicable  to  full-scale  speeds. 

The  next  case  to  be  considered  is  roughness  on  a towed  flat  plate, 

% 

a 21-foot  friction  plane  I}.  The  roughness  was  cast  from  the  bottom  plating 
of  a ship  hull  which  had  been  grit  blasted  to  bare  metal  and  painted  with 
vinyl  paint.  It  was  designated  as  Plate  C in  Reference  2.  It  was  very 
rough  with  a roughness  height  of  about  0.03  inch.  The  drag 
data  is  plotted  in  Figure  4 as.J^  against  log  I^he  A B-characteri zation 

is  calculated  from  Equation  (40)  from  values  of  the  difference  in 
between  the  rough  and  smooth  plates  at  the  same  value  of  Th® 

corresponding  values  of  k*  are  determined  by  using  Equations  (41)  and  (42). 

The  results  are  shown  in  Figure  5 for  theAB-k*  correlation. 

Although  the  flat-plate  roughness  diagram  can  be  calculated  from  the 
A B-characterization,  a more  direct  extrapolation  method  for  the  towed  flat- 
plate  drag  data  will  be  used  to  illustrate  the  method  of  preparing  friction 
lines  for  different  lengths  of  flat  plate.  The  locus  of  a constant AB  is 
shown  in  Figure  6 as  a constant  offset  of  the  smooth  friction  line  in 
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accordance  with  Equation  (37).  A reference  line  for  constant  L*  is 
then  shown  satisfying  Equation  (68)  by  selecting  a convenient  value  of 
L*.  The  loci  of  lines  of  constant  k*  are  shown  offset  from  the  reference 
line  by  a constar.t  distance  in  accordance  with  Equation  (69)  . The 
procedure  is  repeated  for  selected  values  of  the  measured  friction  line. 

Two  lengths  of  300  feet  and  600  feet  are  shown  as  an  illustration.  The 

g 

maximum  value  of  for  the  300-foot  plate  is  7.2  x 10  which  corresponds 
to  a speed  of  about  17  knots.  For  higher  full-scale  speeds,  it  is  then 
necessary  to  increase  the  speed  of  the  towed  friction  plane  to  provide 

extrapolations  for  operating  ship  speeds.  i 

CONCLUDING  REMARKS 

The  three  methods  described  here  of  characterizing  the  drag  increase 
due  to  roughness  have  various  advantages  and  disadvantages.  Hydrodynamically , 
fully-developed  turbulent  pipe  flow  provides  constant  local  conditions  in 
contrast  to  the  variable  turbulent  boundary  layers  on  flat  plates  and 

rotating  disks.  The  three-dimensional  boundary  layers  on  rotating  disks  i 

are  more  complicated  than  the  two-dimensional  boundary  layers  on  flat  plates.  i 

i 

The  difficulty  of  applying  irregular  roughnesses  to  the  inside  of  j 

pipes  may  require  pipe  test  sections  split  longitudinally  for  easy  access.  j 

In  addition,  pipes  have  to  be  long  enough  to  attain  fully  developed  pipe  i 

flow.  Since  the  wall  shear  stress  results  in  a very  small  longitudinal 
pressure  gradient,  the  pipes  have  to  be  very  straight  and  very  level  for 

accurate  measurements.  If  the  roughness  applied  to  the  inside  of  a pipe  ^ 

has  appreciable  thickness,  the  new  effective  inside  diameter  may  have  to  be  1 

accurately  determined. 

A major  problem  with  towed  flat  plates  in  basins  is  obtaining  a 
sufficiently  high  enough  speed  without  undue  free-surface  effects.  However 
this  is  not  a problem  with  plates  in  a water  tunnel.  On  the  other  hand  high 

peripheral  speeds  are  easily  attained  on  a rotating  disk.  In  addition,  a I 

rotating-disk  apparatus  is  usually  a one-purpose  equipment  in  constrast  to  , 

multi-purpose  towing  carriages  or  water  tunnels  used  for  flat  plates  which  j 

compete  for  time  with  other  hydrodynamic  testing.  would  oe  1 

reassuring  to  conduct  definitive  experiments  to  verify  that  the  same  j 

-characterizations  result  irrespective  of  method. 
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APPENDIX  - Enclosed  Rotating  Disk 
1 7 

For  an  enclosed  rotating  disk  the  effective  angular  velocity  is 
reduced  to  ijiu  where  4 is  a factor  accounting  for  the  swirl  of  the 
enclosed  fluid^  0<<(i<l . and  are  modified  accordingly  in  Equations 
(44)  and  (45)  for  use  in  Equations  (46)  through  (51)  as  follows 


and 


A M 


(70) 


(71) 


The  factor  may  be  experimentally  determined  by  a comparison 
of  the  smooth  lines  of  ^ plotted  against  log  R„Vc~  for  enclosed  and 

^ fCjn  ^ ^ 

unenclosed  rotating  disks.  Since  does  not  affect  the  product  R -yc  , 

R ' m 

then  from  definitions 


where  subscripts  en  and  » refer  to  enclosed  and  unenclosed  conditions 

respectively.  The  factor  (p  maybefound^^  to  be  a function  of  R Vc~ 

R m 


(72) 


17  ; — 

Daily,  J.W.  and  R,E.  Nece,  "Chamber  Dimension  Effects  on  Induced  Flow 
and  Frictional  Resistance  of  Enclosed  Rotating  Disks,"  Journal  of 

Basic  Engineering  (Transactions  of  ASNE),  Vol . 82,  Series  D.  pp,  217-252 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  OTIMSRDC  REPORTS,  A FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH- 
NICAL VALUE.  THEY  CARRY  A CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A PRELIM- 
INARY. TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A CASE-BY  CASE 
BASIS. 


